Restoring integrability in one-dimensional quantum gases by two-particle correlations 
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We show that thermahzation and the breakdown of integrability in the one dimensional Lieb- 
Liniger model caused by local three-body elastic interactions is suppressed by pairwise quantum 
correlations when approaching the strongly correlated regime. If the relative momentum k is small 
compared to the two-body coupling constant c the three-particle scattering state is suppressed by 
a factor of (fc/c)^^. This demonstrates that in one dimensional quantum systems it is not the 
freeze-out of two body collisions but the strong quantum correlations which ensures integrability. 

PACS numbers: 67.85.-d,05.30. Jp 



Integrable systems [1] allow for deep conceptual insight 
into many problems of field theory and statistical physics, 
by providing exactly solvable models with rich properties. 
Thermahzation is precluded in integrable models because 
the number of integrals of motion equals to the number of 
degrees of freedom and the 'memory ' of the initial state 
persists in the system forever. A strictly one-dimensional 
(ID) system with local (delta- functional) interactions is 
a prime example of an integrable model [2] . 

In real world, ID systems are realized by strong trans- 
verse confinement {hoj± being the level spacing for the 
radial confinement Hamiltonian), which 'freezes out' the 
radial degrees of freedom. An example are ultracold 
atoms in a tight waveguide where the interval between 
its ground and first radially excited states exceeds both 
the chemical potential {fi <C hu!±) and the temperature of 
the atomic system {ksT <C TiLJi_). Such a system can be 
described by the Lieb-Liniger model jSj. Recently it has 
been shown in experiment that one can cool in ID by 
evaporation to < ■^Tiuj± f3] and theoretically that 

the freeze-out of two body collisions in proportional to 
exp[—21iL0j_/{kBT)] does not necessarily lead to a stop 
of thermahzation [1]. Local three-body elastic collisions, 
where higher radial modes are excited virtually, are only 
suppressed by a ratio of the mean interaction energy to 
?ill!± and lead to thermahzation and a breakdown of in- 
tegrability 0. 

Nevertheless the beautiful experiments by Kinoshita 
et al. [5] showed that for strongly correlated ID systems 
thermahzation is negligible. In this Letter we analyze 
how these two experiments and the theoretical model can 
be reconciled and show that the virtual three-body colli- 
sions described in [3] are suppressed by quantum corre- 
lations in strongly correlated ID Bose gas. 

We start by considering N identical bosons in ID con- 
figuration with the Hamiltonian 

j = l 3 

E U3h{Xj - X.y, Xj - Xj»), (1) 
3>j'>3" 

thereby we use units where Planck's constant is 1 and 



the mass of the atom is | . In these units the strength 
of interaction of two atoms in a tight waveguide is c = 
2as/l\, where Og is the three-dimensional (3D) s-wave 
scattering length and l± ^ is the size of the transverse 
motion ground state of an atom in the waveguide [7] . 

Eq. ([T]) differs from the Hamiltonian of the Lieb- 
Liniger model by including the term U^h that describes 
three-body collisions (its value is non-zero, only if both 
its arguments are approximately 0, i.e., three atoms come 
close to each other), t/ab is obtained by adiabatic elim- 
ination of transverse modes virtually excited by the 3D 
short-range pairwise atomic interaction f?. Note, that ef- 
fective three-body interactions between polar molecules 
emerge in a similar way, due to virtual transitions to an 
off- resonant internal state [8] . 

In a non-degenerate state, where free phase space is 
available for the scattered particles, J/sb leads to ther- 
mahzation due to scattering of particles, since a three- 
particle elastic collision in ID, unlike a two-particle col- 
lision, allows for redistribution of the colliding particles 
momenta. In the ground state, this type of interactions 
not only changes the energy of the system, but also pre- 
cludes the Lieb-Liniger quasi-momenta from being inte- 
grals of motion. U^^ is the source of thermahzation and 
leads to a beakdown of integrability. 

In what follows, we calculate this three-body scattering 
amplitude in the presence of the delta-functional pairwise 
interactions. The real part of this scattering amplitude 
determines the vertex for the effective three-body inter- 
action by analogy with two-particle scattering in two di- 
mensions m Uni E] . The thermahzation rate IS propor- 
tional to the square of its absolute value. The stronger 
the pairwise interparticle repulsion, the smaller is the 
probability of a close encounter of three particles. This 
results in suppression of the three-body scattering ampli- 
tude, which is the main subject of the present Letter. 

To analyze the Hamiltonian ([T|) for iV = 3 we express 
it in hyperspherical coordinates R, a defined as |12j 



_Rsin a 



Xi - X2 



V2 



R cos a 



X3 



Xi + X2 

2 



(2) 



2 




FIG. 1: Transformations of the hyperangle a corresponding 
to the pairwise permutations of particles j <-> j' . The regular 
hexagon indicates the symmetry of the Hamiltonian of three 
pairwise-interacting identical particles. The two-body inter- 
action potential has delta-functional singularity on straight 
lines (not shown in this figure) going through opposite ver- 
tices of the hexagon, cf. Eq. ([sjl. See text for the details. 



and the center-of-mass coordinate X = \{xi + X2 + x^): 
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If one deals with three particles in the 3D space, the 
total number of degrees of freedom equals to 9. After 
separation of the center-of-mass motion, six degrees of 
freedom remain. The corresponding coordinates in the 
usual 3D case are R, a and four angles characterizing the 
orientation of two independent differences of the particles 
radius- vectors [T5]. In our ID case, the co-ordinate space 
reduces, after separation of the center-of mass- motion, to 
a plane, R and a being its polar co-ordinates. 

For indistinguishable bosons the permutation xi <-> X2 
corresponds to the reflection operation that changes a to 
~a. The two remaining pairwise permutations, X3 <-*■ X2 
and Xi 2:3, correspond to the product of this reflection 
and rotation in the plane by — 27r/3 and 27r/3, respec- 
tively, as shown in Fig. 1. The two cyclic permutations 
of xi, X2, X3 correspond to mere rotations by ±27r/3. 
The three-body interaction potential (its particular form 
will be specifled later) must be invariant with respect to 
the permutation of the particles. The Hamiltonian of the 
three particles interacting via pairwise interactions only, 
i.e., Eq. ([3| except of the term [/3b(i?, a), possesses a 
richer symmetry (the symmetry group Dg of a regular 
hexagon [T3] ; two of its six reflection-symmetry axes are 
shown in Fig. 1 by dash-dotted lines). Adding any rea- 
sonable effective three-body interaction does not change 
this symmery. 

The Schrodinger equation for the three-particle wave 
function is H'i>{xi,X2,X3) = {kf -I- -I- k^)'if{xi,X2, X3). 
The wavenumbers kj can be defined from the set of 
transcendental equations [2], provided that the periodic 



boundary conditions are set on the interval of the length 
L. By setting L — > cx), we obtain a continuous spec- 
trum of real kj's in the repulsive case (c > 0). Then we 
separate the center-of-mass motion and describe the rela- 
tive motion in hyperspheric coordinates: ^'(xi, a;2, xa) = 
exp[i{ki + k2 + A:3)X]?/'r(-R, «)• The kinetic energy of the 
relative motion is k^ = |[(A:i — fc2)2 -|- {k2 — k^Y + {k^ — 
kiY]. We decompose the wave function of the three col- 
liding particles in the center-of-mass frame as |12j 



V'r(-R,a) = ^F„(i?)x„(i?,«), 



(4) 



n=0 



where Xn{R,ot) is an eigenfunction of the hyperangle- 
dependent part of the Hamiltonian: 



V2cR J2 ^( 



\l{R)Xn{R.a). (5) 



i/' = -1,0,1. (6) 



Due to the bosonic symmetry, 

Xn{R, a) = XniR, ±a + 2nv'/3), 

Eq. ^ defines the spectrum of eigenvalues and, hence, 
the adiabatic hyperspheric potentials \\{R)/R^. The 
eigenfunctions are 

X„(i?,a) = Xn(^, a)+Xn(^, a-27r/3)+x„(i?, a+27r/3), 

(7) 

where 



Xn{R, a) 



cos[A„(i?)(7r/6 - |a|) - 7rn/2], |a| < 7r/3 
0, otherwise 

(8) 

and the eigenvalues are the positive roots of the equation 



A„(i?)tan[7rA„(i?)/6-7rn/2] = ci?/V2, n = 0, 1,2, .... 

(9) 

After integrating out the hyperangular variable, the 
Schrodinger equation reduces to the set of coupled ordi- 
nary differential equations 
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F„ = fc2F„,(10) 
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Neglecting all ofF-diagonal (coupling) terms is termed 
adiabatic hyperspherical approximation |14j . However, 
for the sake of clarity, we neglect also diagonal terms 
Wnn and Ynn [12 . Indeed, the respective terms decrease 
more rapidly than others by a factor of as i? — > cxd, 
and at ci? < 1 the correction to the adiabatic potential 
\q{R)/R^ due to Yqq{R) is less than 5 %. The finite value 
of Woo can be neglected compared to Xq{R)/R^ — > oo as 
0. 

In what follows we assume the long-wavelength limit 
k <^ c, which holds in the regime of strong interactions [2] 
at temperatures of the order of or less than the chemical 
potential. It enables us to consider in Eqs. (|4] [lO| only 
the lowest hyperangular mode, n — 0, because in higher 
modes scattering is not efficient. We introduce the effec- 
tive radius tq of the three-body interactions: U^h rapidly 
vanishes for interatomic separation values exceeding vq. 
For atoms trapped in a tight waveguide, tq, like the ef- 
fective range of two-body collisions, is much smaller than 
other available length scales, in particular, vq ^ c^^. We 
consider now, under the assumptions mentioned above, 
the equation 



-—R—F 
RdR dR ° 



i?2 



UooiR) 



Fo = A;2Fo (15) 



with the boundary conditions requiring Fq to be finite for 
both R — and i? — > oo, for the "partial wave" corre- 
sponding to the lowest eigenvalue Ao(i?), whose asymp- 
totic expressions are 
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(16) 



We assume that the averaged over the hyperangle three- 
body interaction potential Uqo{R) vanishes at R > vq, 

where tq <C c~^. 

First we solve Eq. (15) for R > tq. For ci? < 1 the 



two independent solutions are the modified Bessel func- 
tions of the zeroth order [T5] Iq{2^/PR) and Ko{2y/PR), 
where /3 = (3'\/2/7r)c (we call them the inner solutions). 
The outer real solutions hold for cR ^ 1 and are the 
Bessel functions of the third order J3{kR) and Y3{kR). 
Now we note, that in the range ^ R <C k~^ the 
outer solutions can be represented by their asymptotics 
J3{kR) « {kRf/48 and Y^ikR) « (-16/7r)/(fci?)3. The 
outer and inner solutions can be tailored by quasiclassical 
solutions in two dimensions 



/'qc(^) 



exp 



± 



^R^'-fl 



(17) 



Note that two-dimensional quasiclassical solutions differ 
from one-dimensional ones [16^ by a prefactor 1/\/R. By 
a proper choice of the constants we match at i? ~ 
the quasiclassical solutions to the inner solutions: <f)^^{R) 
to Io{2y/jm) « (47rV^)-i/2exp(2V;3i?) and 0qc(i?) to 
Ko{2^/PR) « (4V;gi?/7r)-i/2exp(-2V;5i?), respectively. 



On the other hand, the asymptotics of the quasiclassical 
solutions at ci? ^ 1 are 



<c(^) 



(i?) 



12S 
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where S w 0.18 is a numerical constant. This tailoring 
enables us to construct two independent real solutions 
F^{R) of Eq. (15) in the whole range R > tq with the 
following asymptotics: 

j WVM), cR<l 
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(20) 

where & is a real number of the order of or less than 1 . In 
other words, the term Y^lkR) oc R~^ is only the leading 
term of the expansion of F^ (R) in the intermediate range 
of R; the solution J2,{kR) vanishing for small kR <C 1 is 
contained, in a general case, in terms beyond the accu- 
racy of the quasiclassical approximation |16] . The upper 
limit to the magnitude of h follows from the fact that 
the quasiclassical expression (18) for (t)~^{R) holds up to 
k-\ _ 

for i? > ro is then 



R 



The general solution of Eq. ( 15 



Fo{R) = C^F+{R) + C2F-{R). 



(21) 



The particular ratio between the coefficients Ci and C2 
is found from matching the logarithmic derivatives of the 
solution given by Eq. (21 1 and its solution at R < rg. 

As a concrete example, illustrating the basic physics 
through simple analytic expressions, we take the three- 
body interaction potential in the form of the sum of a 
rectangular potential well and a term compensating the 
adiabatic hyperspherical potential at _R < tq: 



UooiR) 



P/R, 



R < ro. 



(22) 



Without loss of generality, we assume further that k ^ q. 
In this case the regular solution at i? < tq is Jo{qR). 
Matching the logarithmic derivatives of the solutions at 
R > ro and R < ro, recalling that ro <C , and using 
the asymptotics Io{z) « 1 and Ko{z) « — Inz at 2: ^ 0, 
we obtain 



C2 



= ln(2^ 
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(23) 



Then we calculate the partial scattering amplitude /o 
(corresponding to the scattering channel with n = 0) 
introduced by expressing Eq. (21 1 at i? ^ cx) in the form 



Fo{R) w const [J^ikR) 



i}oH\^\kR)], 



(24) 



where Hi^ '{z) = J3{z) + iY^^z) is the Hankel function of 
the first kind, corresponding to the outgoing (scattered) 
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wave. Note that there are various definition of the scat- 
tering amplitude in 2D, diS^ering by a complex prefactor 
[TUl [Tfl [15] . The rate of three-body collisions is propor- 
tional to I /o lo- 
using Eqs. (19 20), we obtain 



/o = s ^- (25) 

where fl = 384(S/7r)2 w 1.26 and J = 

>^o('/''o)/[9''o>^i(9''o)]- Other choice of f/oo(-R) rather than 
Eq. ( p2| gives different value for J^, but does not change 
the general structure of Eq. (251. Although the exact 



value of b can not be obtained from the quasiclassical 
matching of solutions, it is relatively small (|6| < 1) and 
thus b can be neglected in Eq. ( 25 ) . 

From Eq. ( 25 1 we conclude that the three-body scat- 
tering amplitude decreases in proportion to (fc/c)^ as 
c — *■ cxD, and the three-body scattering rate in a ID 
system of bosons in the case of strong pairwise inter- 
action is suppressed by a factor ^ (fc/c)"'^^. This suppres- 
sion is the main result of the present paper. Averaging 
over collision momenta in a moderately-excited strongly- 
interacting state we obtain the scattering rate suppres- 
sion factor ~ 7~^^, where 7 = c/nin is the Lieb-Liniger 
parameter [5] and ni£i = {'^^ {x)'^ {x)) is the ID density 
of particles, ^(x) being the bosonic field annihilation op- 
erator. 

We can now compare our result with the zero- 



distance three-particle correlation function 173(0) = 
{'^^^{x)'^^{x)) /n\j-,. In the strong interation limit 7 ^ 
1, where 53(0) cx 7"^ [13 [20]. This thus proves the con- 
jecture [1] that the pairwise interactions and the quantum 
correlations induced by them in a strongly-interacting ID 
bosonic system suppress the three-body elastic scatter- 
ing rate, and, hence, thermalization, by a factor oc 53(0). 
In other words, strong quantum correlations restore inte- 
grability as a ID system approaches the Tonks-Girardeau 
regime. This observation is in agreement with the exper- 
imental results [6; . 

In conclusion we remark, that integrability in tightly 
confined ID systems is not, as suggested by quasiclas- 
sical arguments, caused by the freeze-out of two-body 
collisions, but is ensured by the quantum correlations in 
a strongly interacting ID system. The genuinely quan- 
tum effect of virtual excitations to higher radial states 
opens a way to thermalization by three-body collisions 
and breaks integrability in the quantum description of 
a ID system. These three-body collisions can be sup- 
pressed by quantum correlations caused by strong pair- 
wise repulsions. If they dominate, as in a strongly cor- 
related ID Tonks-Girardeau gas, they restore the inte- 
grability of the system. This points to an interesting 
difference between classical and quantum systems. 
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